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The phase diagram of the massive Gross-Neveu model, revisited 
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The massive Gross-Neveu model is solved in the large A'' limit at finite temperature and chemical 
potential. The phase diagram features a kink-antikink crystal phase which was missed in previ- 
ous works. Translated into the framework of condensed matter physics our results generalize the 
bipolaron lattice in non-degenerate conducting polymers to finite temperature. 
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In its original form, the Gross-Neveu model jll is a rela- 
tivistic, renormalizable quantum field theory of A'^ species 
of self-interacting fermions in 1+1 dimensions with La- 
grangian 
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i=l 
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The bare mass term ~ mo explicitly breaks the discrete 
chiral symmetry ■0 — > ^^ip of the massless model. As 
far as the phase diagram is concerned, the 't Hooft limit 
TV^cxD, (7^~l/iVis most interesting since it allows one 
to bypass some of the limitations of low dimensions, while 
justifying a semi-classical approach. Interest from the 
particle physics side stems from the fact that the simple 
Lagrangian shares non-trivial properties with strong 
interaction physics, notably asymptotic freedom, dimen- 
sional transmutation, meson and baryon bound states, 
chiral symmetry breaking in the vacuum as well as its 
restoration at high temperature and density (for a re- 
view, see e.g. 0)- It also has played some role as a 
testing ground for fermion algorithms on the lattice 
Perhaps even more surprising and less widely appreciated 
is the fact that Gross-Neveu type models have enjoyed 
considerable success in describing a variety of quasi-one- 
dimensional condensed matter systems, ranging from the 
Peierls-Frohlich model Q over ferromagnetic supercon- 
ductors [3 to conducting polymers, e.g. doped trans- 
polyacetylene By way of example, the kink and kink- 
antikink baryons first derived in field theory in the chiral 
limit (mo = 0) have been important for understand- 
ing the role of solitons and polarons in electrical con- 
ductivity properties of doped polymers . Likewise, one 
can show that baryons in the massive Gross-Neveu model 
(mo ^ 0) are closely related to polarons and bipolarons 
in polymers with non-degenerate ground states, e.g. cis- 
polyacetylene ^8, ,9]. Finally, Gross-Neveu models with 
finite N have been found useful for describing electrons 
in carbon nanotubes and fuUerenes jlO|. 

The claim that a relativistic theory is relevant for con- 
densed matter physics is at first sight very provocative 
|ll| . A closer inspection shows that a continuum approx- 
imation to a discrete system, a (nearly) half-filled band 
and a linearized dispersion relation of the electrons at 
the Fermi surface are the crucial ingredients leading to a 
Dirac-type theory. The Fermi velocity plays the role of 



the velocity of light and the band width the role of the UV 
cutoff. In some cases the correspondence is so close that 
results can be taken over from one field into the other. 
This was illustrated in where we borrowed results 
from the theory of non-degenerate conducting polymers, 
in particular the bipolaron lattice, for solving the zero 
temperature limit of the massive Gross-Neveu model. 

The Lagrangian ^ has two bare parameters, and 
mo. In the process of regularization and renormalization 
all observables can be expressed in terms of two physical 
parameters m and 7. The relation to the bare quantities 
and the UV cutoff A is given by the vacuum gap equation 
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Whereas m merely provides the overall mass scale and 
can be set equal to 1, the renormalized fermion mass ra- 
tio 7 (called "confinement parameter" in condensed mat- 
ter physics) parametrizes different physical theories. It 
measures the amount of explicit chiral symmetry break- 
ing and vanishes in the massless (mo — 0) case. The 
phase diagram of the model (QJ was proposed in 1995 
[l^. but recent findings showed that the result is flawed. 
A related problem in the massless Gross-Neveu model 
was pointed out and subsequently cured by us in pre- 
vious works [13L IT^ ITHj l . We will recover these findings 
below as special case 7 = of the present work. 

Let us first recall the original analysis !l2!| which is still 
valid in a certain region of the phase diagram. In the 
large A'' limit, a saddle point approximation to the func- 
tional integral or, equivalently, a (thermal) Hartree-Fock 
approach becomes exact. One finds that the fermions 
condense to form a scalar potential S{x). The authors 
proceed along similar lines as the 1985 work at 7 = 
and tacitly assume that S is spatially constant lead- 
ing to a (/i, T)-dependent dynamical fermion mass M . 
This effective mass is determined via self-consistency by 
the thermal expectation value of the fermion condensate. 



M = mo-A^g2(i^V)th. 
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For a translationally invariant system, this in turn is tan- 
tamount to minimizing the renormalized grand canonical 
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(to = !,£■ = \/ q'^ + Af 2) with respect to Af. Depend- 
ing on the parameters {fi, T, 7), there may be one or two 
local minima with the possibility of a first order phase 
transition. The phase diagram depending on /i, T, 7 is 
plotted in Fig. 1. In the chiral limit (7 = 0), the mas- 




FIG. 1: Phase diagram of the massive Gross-Neveu model as 
a function of 7, /i, T, (incorrectly) assuming unbroken transla- 
tional invariance. The various phase boundaries are explained 
in the main text. 

sive phase at low (/i, T) is separated from the chirally 
restored, massless phase at high (/i, T) by a critical line 
(APtB) The upper part of this line (APt) is 2nd 

order, the lower part (PtB) first order, with a tricriti- 
cal point Pt separating the two. As the parameter 7 is 
switched on, the 2nd order line disappears in favour of a 
cross-over where the fermion mass changes rapidly, but 
smoothly. The first order line on the other hand survives, 
ending at a critical point. If plotted against 7, these crit- 
ical points lie on the third curve (PtC) emanating from 
the tricritical point shown in Fig. 1. For 7 > the effec- 
tive fermion mass M is different from zero everywhere. 
If one crosses the shaded critical "sheet" in Fig. 1, the 
mass changes discontinuously, dropping with increasing 
chemical potential. More details on the thermodynamic 
implications of such a phase diagram are given in 12]. 

In the meantime, evidence has accumulated that this 
picture cannot be correct. This has so far only been 
clarified in the two coordinate planes 7 = and T = 
^3] of Fig. 1, and we propose to complete the phase 
diagram here. Briefly, what has been missing is a proper 
treatment of the baryons which can give rise to a kink- 
antikink crystal phase in certain regions of temperature 
and chemical potential. Actually, this phenomenon has 
been known in condensed matter physics for a long time: 
It reflects the Peierls instability, a general phenomenon 
in one-dimensional systems. In order to lower the energy 



or the thermodynamic potential, a periodic modulation 
of the lattice is induced, thereby opening a gap at the 
Fermi surface. In the relativistic Gross-Neveu model it 
turns out that the scalar condensate {"ipip) (the a-field) 
plays the role of the ion displacement. 

To incorporate this kind of physics into the Gross- 
Neveu model, one has to solve the Dirac equation self- 
consistently (including the negative energy sea) with a 
spatially varying, periodic scalar potential ^(a;), in gen- 
eral a highly non-trivial task even in one space dimen- 
sion. Fortunately, in the special cases considered so far 
(namely 7 = and T = 0), this has turned out to be ex- 
actly feasible due to a number of lucky circumstances. 
The self-consistent potentials are such that the Dirac 
equation can be mapped onto a 2nd order Schrodinger 
type equation with the Lame potential of order 1. The 
single particle eigenspinors are known in terms of theta 
functions. This enables one to perform an exact calcula- 
tion of both the ground state and the phase diagram, in- 
cluding an analytical proof of self-consistency. At 7 = 0, 
such a calculation was done independently in condensed 
matter physics P| and relativistic field theory [l^ . Dur- 
ing the last two decades, interest in non-degenerate con- 
ducting polymers has fostered a lot of theoretical work on 
the polymer side. As a result, analytical bipolaron lat- 
tice solutions at zero temperature were found by several 
groups 

CI III USUI El- We have recently shown that 
the x-dependent gap parameter of the bipolaron lattice 
is adequate for solving the T = 0, massive Gross-Neveu 
model fl^ . 

We have extended this calculation to finite tem- 
perature and again found exact self-consistency. We 
can therefore now present a revised phase diagram in 
the 3-dimensional (fi, T, 7)-space which supersedes the 
old, translationally invariant phase diagram of Fig. 1. 
Whereas the solution demands some mathematics (simi- 
lar but more complex than the 7 = solution presented 
in C3) the outcome can be visualized in a single graph. 
This made us decide to present in this Letter our main re- 
sults without derivation. The detailed formalism and fur- 
ther results on the revised thermodynamics of the mas- 
sive Gross-Neveu model will be given elsewhere. 

The crucial observation is that the scalar potential at 
T ^ has the same general functional form as at T = 0. 
It can be written in terms of three parameters A, b, k and 
Jacobi elliptic functions of the modulus k (suppressed 
here) as 
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This form equals the more complicated expression em- 
ployed in the literature on the bipolaron lattice and in 
our T = paper C3- As mentioned above, the Dirac 
equation with potential © is equivalent to the single gap 
Lame equation. Using the known spectrum and density 
of states of the Lame potential, one can construct the 
grand canonical potential and minimize it with respect 
to {A, b, k). Due to the specific structure of the equations 
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obtained, it is possible to reduce the solution to basically 
solving a transcendental equation in one variable contain- 
ing a one-dimensional numerical integration. The proof 
of the self-consistency condition which replaces Eq. , 

S{x)^mo-Ng^{'4,i;)ti,, (6) 

establishes that the minimum of the grand canonical po- 



tential is an exact solution of the large TV Gross-Neveu 
model. Comparison with the translationally invariant 
solution shows that the crystal is thermodynamically 
favourable whenever it exists. 

We now turn to the results, focussing on the phase 
boundaries as depicted in Fig. 2. 
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FIG. 2: Revised phase diagram of the massive Gross-Neveu 
model. The shaded surfaces 1, II separate the kink-antikink 
crystal from a massive Fermi gas and correspond to 2nd order 
phase transitions. 



In the 7 = plane, the result for the revised phase di- 
agram including the kink-antikink crystal is known (15|: 
The 2nd order line (APt) of the old phase diagram is un- 
affected. The first order line (PtB) of Fig. 1 is replaced 
by two 2nd order lines delimiting the kink-antikink crys- 
tal phase. The tricritical point is turned into another 
kind of multicritical point at the same (/x, T) values. As 
we turn on 7, the second order line separating massive 
[M > 0) and massless (M — 0) phases disappears as 
a consequence of the explicit breaking of chiral symme- 
try. The crystal phase survives at all values of 7, but 
is confined to decreasing temperatures with increasing 7. 
For fixed 7, it is bounded by two 2nd order lines join- 
ing in a cusp. The cusp coincides with the critical point 
of the old phase diagram but has a significantly different 
character. The crystal phase exists and is thermodynam- 
ically stable inside the tent-like structure formed out of 
two sheets denoted as / and //. These sheets are de- 



fined by K = (/) and k = 1 {II), respectively. The 
line where they join corresponds to & = and coincides 
with line (PtC) in Fig. 1. The baseline of sheet // in the 
(/i,7)-plane has a simple physical interpretation: It re- 
flects the 7-dependence of the baryon mass in the massive 
Gross-Neveu model, or equivalently the critical chemical 
potential at T = 0. The chiral limit 7^0 can be iden- 
tified with b K(k). In this way, one can recover the 
simpler 2-parameter ansatz for the self-consistent poten- 
tial used in 15] from the 3-parameter ansatz Eq. (jsj. 

In order to further illustrate the nature of the phase 
transitions between inhomogeneous and homogeneous 
phases, let us remark that outside the tent the scalar 
potential is ^-independent, i.e. a dynamical mass. If one 
crosses the tent from side I to side // on some curve, 
the parameter k varies continuously from to 1. By way 
of example consider the straight line 7 = 0.5, T — G.05/i 
which pierces sheet I ai ^ — 1.1 and sheet II at 11 = 0.85. 
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FIG. 3; Evolution of self-consistent S{x) as one crosses the 
crystal region on a straight path. As fi decreases from 1.1 to 
0.85, the elliptic modulus k varies between and 1. 



Just outside of the crystal phase, the dynamical mass val- 
ues we find are 0.4 (/) and 1.0 (//). According to the old 
phase diagram, the mass jumps suddenly from 0.54 to 1.0 
upon crossing the first order sheet at = 0.87. By con- 
trast. Fig. 3 shows the continuous evolution of S{x) in the 



exact calculation. At k = 0, an instability occurs with 
respect to oscillations of a finite wavenumber. As k in- 
creases, the amplitude of these oscillations become larger 
whereas the period first grows rather modestly. When k 
approaches the value 1 the period grows rapidly and, save 
for small dents representing widely spaced baryons, S{x) 
reaches a constant value connecting to the translation- 
ally invariant phase. At k = 1 the system is instable 
against single baryon formation. This subtle interpola- 
tion of S{x) between two constants caused by the Peierls 
instability is only crudely modeled by the translationally 
invariant scenario. 

Concluding, we find it gratifying that the simple 
Lagrangian gives rise to such a rich phase diagram. 
We emphasize that we would not have been able to solve 
this problem without the works on the bipolaron lattice 
in conducting polymers. We hope that our extension 
to finite temperature will in turn lead to new results in 
condensed matter systems. 

We should like to thank Wolfgang Hausler for helpful 
conversations on one-dimensional condensed matter 
physics. 
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